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We study the response to shear of a one-dimensional monolayer embedded between two rigid 
' plates, where the upper one is externally driven. The shear is shown to excite "shearons", which 

are collective modes of the embedded system with well defined spatial and temporal pattern, and 
which dominate the frictional properties of the driven system. We demonstrate that static friction, 
stick-slip events, and memory effects are well described in terms of the creation and/or annihilation 
of a shearon. This raises the possibility to control friction by modifying the excited shearon, which 
■ we examplify by introducing a defect at the bottom plate. 

m 

i 1 PACS numbers: 46.55.+d, 81.40.Pq, 68.15.+e, 05.45. -a 

^ \ The field of nanotribology evolves around the attempts to understand the relationship between macroscopical 
frictional response and microscopic properties of sheared systems |l|] . New experimental tools such as the surface force 
apparatus (SFA) are used to explore shear forces between atomically smooth solid surfaces separated by a nanoscale 
molecular film [E[. These experiments have unraveled a broad range of phenomena and new behaviors which help 
shed light on some "old" concepts which have been considered already textbook materials: static and kinetic friction 
forces, transition to sliding, thinning, and memory effects. These and other observations have motivated theoretical 
efforts [p|3,P|,P|,[^|,PP JlO| , pd| , P^ | , both numerical and analytical, but many issues have remained unresolved, in particular 
the relation between the macroscopic observables and the microscopic properties of the embedded system. 

In this Letter we introduce the concept of "shearons", which are shear-induced collective modes excited in the 
| embedded system and characterized by their wave vector q. Shearons, which display well defined spatial and temporal 
c*j . patterns, dominate the frictional properties of the driven system and are found useful in establishing a connection 
between frictional response and motional modes of the embedded system. Within this framework observations such 
as static friction, stick-slip behavior, and memory effects can be correlated with the creation and/or annihilation of 
a shearon. These correlations suggests the possibility to control friction by modifying the shearon's wave vector and 
t— ( \ thereby tuning the embedded system, by adding, for example, a defect at one of the plates. 

We start from a microscopic model which has been investigated recently and has been shown to capture many of 
the important experimental findings |||^|. The model system consists of two rigid plates, with a monolayer of A" 
particles with masses m and coordinates Xi embedded between them. The top plate with mass M and center of mass 
coordinate X is pulled with a linear spring of spring constant K. The spring is connected to a stage which moves 
, with velocity V. This system is described by A^ + 1 equations of motion 

N N „ , _ „, 

MX + ^v(X- x t ) + K(X -Vt) + ^2 g X = (1) 
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The second term in Eqs. (|l|) and (||) describes the dissipative forces between the particles and the plates and is 
proportional to their relative velocities with proportionality constant r/, accounting for dissipation that arises from 
interactions with phonons and/or other excitations. The interaction between the particles and the plates is represented 
by the periodic potential Q(x) = — $o cos(2-7ra/6). Concerning the inter-particle interaction, we assume nearest 
neighbor interactions of two types: (i) harmonic interaction ^(xj — Xi±%) = (k /2) [xj — Xi±i ± a] 2 and (ii) Lennard- 
Jones interaction ^{xi — Xi±\) — (ka 2 /72){[a/(xi — Xi±\)] 12 — 2\a/(xi — Xi±i)] 6 } p3| . The two plates do not interact 

directly. 

The basic frequency is chosen to be the frequency of the top plate oscillation in the potential fl = (2w/b)y/ N&q/M . 
The other frequencies in the model are the frequency of the particle oscillation in the potential ui = (2tt/&) -\/<&o/m, 
the characteristic frequency of the inter-particle interaction lo = yjk/m, and the frequency of the free oscillation of 
the top plate ft = \J K/M. To simplify the discussion we introduce unitless coordinates Y = X/b and yi = Xi/b of the 
top plate and the particles, respectively, as well as an unitless time r = fit. We define the following quantities: The 
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misfit between the substrate and inter-particle potentials' periods A = 1 — a/6, the ratio of masses of the particles and 
the top plate e = Nm/M, the unitless dissipation coefficient 7 = Nr]/(M£l), the ratio of frequencies of free oscillation 
of the top plate and the oscillation of the top plate in the potential a = fi/fi, the ratio of the frequencies related to 
the inter-particle and particle/plate interactions f3 = lo/lj, and the dimensionless velocity v = V/(bQ,). Additionally, 
the friction force per particle /k is defined as /k = F^/(NF S ), where is the total friction force measured with the 
external spring and F s = 2n& /b is the force unit given by the plate potential. Here, we keep fixed the number of 
particles N — 15, the mass ratio e = 0.01, the misfit A = 0.1, and the relative strength of the inter-particle interaction 
f3 2 = 1. We vary only the relative strength of the external spring a 2 , the dissipation coefficient 7, and stage velocity 
v. 

In order to analyze the motion of the embedded system more closely, we seperate the motion of the particles into 
the center of mass part y cms = 1/Aj^^ Hi [0 an d the fluctuations Syi = yi — y C ms- It has been observed in similar 
models |l^Jl^,|l^] that different modes of motion can coexist for a given set of parameters and lead to different frictional 
forces. Here, we concentrate on these solutions of the coupled dynamical equations (0) and (0) which correspond to 
smooth or to stick-slip motion of the top plate. To understand the nature of the motion of the embedded system in 
these regimes and the relation to the frictional response, we choose the particles density g as a observable. Instead of 
defining the density as X)i=i S(u ~ Hems ~ Syi), we represent each particle by a Gaussian of width a — 1, namely 



N 



g(y - y cms ,T) = ^exp • 



y - Vcms - Syi 



(3) 



This allows to visualize correlated motions of the particles on length scales which exceed nearest neighbor distance 
and are of the order of 2a to 3a. Using Eq. (||) we find that the response of the embedded system to shear can 
be described in terms of collective modes which result in well defined spatial and temporal pattern in the density, 
which we call "shearons". These shearons are the spatial/temporal manifestation of parametric resonance between 
the external drive and the embedded system |l7|,[l8| . 

In Fig. 0(a)-(c) we present three stable shearons for a chain of harmonically interacting particles with free boundary 
conditions, where in all cases the top plate slides smoothly. All three shearons share the same set of parameters and 
differ only in their initial conditions. We find that the observed mean friction force per particle (/k) ((•} denotes time 
average), which is directly related to the spatial/temporal fluctuations by |Q 

<A)=^ + ^|>A , (4) 

decreases with increasing shearon wave vector g. This means that the mean friction force can be reduced by increasing 
the shearon wave vector. In the present example, we observe a reduction from (/k) = 0.263 [Fig.0(a)] to (/k) — 0.180 
[Fig.0(b)], and to (/k) = 0.125 [Fig.0(c)]. The lowest possible mean friction force (/k) = n^fv can be achieved in the 
limit q — > 00. For our choice of parameters n^/v = 0.04712. 

We find that the regions of high density fluctuations exhibit low velocity fluctuations and, according to Eq. (0), low 
dissipation. These regions are surrounded and stabilized by regions of low density fluctuations with highly fluctuating 
velocities and hence high dissipation. We note that the minima/maxima of the density reflect correlations in the 
motion of the neighboring and next-neighboring particles rather than actual locations. The cooperative nature of 
motion is visualized here introducing a finite width a. 

As we have already seen, larger wave vectors correspond to higher average density of the embedded system. This 
means that for a given number of particles the mean chain length (L) = {y^ — yi) decreases with q: (L) = 13 
[Fig. 1(a)], (£) = 12.5 [Fig. 0(b)], and (L) £S 12 [Fig. 0(c)]. Note that for N = 15 and A = 0.1 the equilibrium length 
of the free chain is 12.6. The interplay between mean length (L) and the shearon wave vector q becomes clearer 
within the calculations under periodic boundary conditions, where q is a function of the box size £, see Fig. 0(d)- 
(f). In Fig. 0(d) we present the spatial/temporal pattern found for harmonically interacting particles within a box 
of size £ = 13. It is found that for similar wave vectors the mean friction force is smaller for the case of periodic 
boundary conditions compared to a free chain [c.f. Fig. 0(c)], since the free ends strongly fluctuate and cause additional 
dissipation. 

When replacing the harmonic with the Lennard- Jones interaction [see Fig. 0(e) and (f)], the picture remains 
basically unchanged. The essential differences are found for: (i) very large box sizes, when the embedded Lennard- 
Jones system breaks and one reobtaines the independend particle scenario described in ||; (ii) very small box sizes, 
when the hardcore of Lennard- Jones potential becomes dominant. The latter case corresponds to the high pressure 
limit p oc i~ l —> 00 and results in a strong stiffness. This, together with incommensurability, allows to achive the 
lower bound for the mean friction force (/k) = tt^/v, c.f. Eq. (0). In the case of harmonically interacting particles a 
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minimum of the mean friction force (/k) as a function of box size £ is found for a finite value of £, which is determined 
by commensurate condition of shearon wave vector and the box size. These effects will be discussed elsewhere |]l9| . 

We now explore the concept of shearons in relationship to varios frictional phenomenon. We start from the well 
known stick-slip phenomenon observed in many nanoscalc systems at low driving velocity, but whose nature is still 
not well understood. The start of a slip event has been commonly attributed to the 'melting' of the embedded system, 
namely a transition from a ordered 'solid'-like to a disordered 'liquid'-like structure, that 'refreezes' at the end of 
the slip event. In Fig. |^(a) we show that during slippage the motion of the embedded system is highly ordered and 
highly correlated. At the start of the slip event, the moment of the highest spring force [enlarged in Fig. §(b)], a 
shearon is created, persisting with a constant wave vector until it gets annihilated at the moment of the lowest spring 
force [enlarged in Fig. §(c)]. The shearon gets annihilated since it cannot exists below a certain shear force needed 
to compensate the energy dissipation. As a result we find that the static force / s needed to be overcome in order to 
initiate the motion is the shear force needed to create the shearon. This becomes clearer in a stop/start experiment, 
where for a smoothly sliding top plate the external drive is stopped for a certain time and reinitiated afterwards. We 
find that the static friction force needed to restart the motion manifests a stepwise behavior as a function of stopping 
time |p|. It vanishes only as long as the motion is restarted within the lifetime of the shearon, giving a possible 
explanation of the memory effects observed in [^0| . 

Since in the shearon description the frictional force is determined by the wave vector, it is suggestive to influence 
the friction by modifying q. It is possible to change the shearon wave vector by changing external parameters such as 
the stage velocity (a higher velocity results in general in a larger wave vector). A shearon with a large wave vector, 
created at a high velocity, can be maintained at lower velocities by deceleration, giving rise to a hysteretic behavior 
observed in many experimental systems . 

Here, we present as an example a "chemical" method to manipulate the shearon wave vector by introducing a 
defect. The defect is placed at the bottom plate at an integer position yo, leading to a modified bottom plate 
potential &y {y) — $0(1 — h{l + cos(27r[j; — Vo})}) for \y — yo\ < 1/2 and ^' yo {y) = $(2/) otherwise. Here h defines 
the relative depth of the minimum at position yo, with h = 1 being the regular case. As an example, the resulting 
density pattern for h = 1/2 is shown in Fig. |l where it can be seen that scattering by the defect changes the shearon 
wave vector q into a new wave vector q' > q. The new shearon with wave vector q' is stable, leading to a decrease of 
the mean friction force from (/k) = 0.263 before passing the defect to (/k)' — 0.178 afterwards. Depending on the 
amplitude of h, both a decrease and an increase of the friction force are possible, which provides a method to tune 
(/k). It was already observed in similar models that disorder can significantly change the frictional behavior ||T^ , |l6|| . 

In this Letter we have introduced the concept of shearons, which are shear-induced collective modes in the embedded 
system, and have demonstrated that their properties dominate the frictional behavior of the driven system. The results 
have been obtained for a certain set of parameters describing the embedded system. However, the collective effects 
discussed above have a general nature and have been found in a wide range of parameters, except for the limits 
N — > 1, /3 — s- 0, and/or A — > 0, which correspond to the case of independent particle motion. We believe that the 
observed collective modes of motion are even more general and should not depend on the one-dimensionality of our 
model. In particular, we expect shearons to persist in systems of higher dimension, as long as the embedded system 
remains a monolayer, so that each particle interacts with both the bottom and the top plate. The possibility to modify 
shearons and hence the frictional response by (i) (ambient) pressure and (ii) defects at the plates should realizable 
experimentally. 
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FIG. 1. Plot of the particles density g vs position y — j/ cms and time r for different shearons. The white and black colors 
indicate low and high density, respectively [the gray scales are chosen independently for each subfigure to maximise contrast]. 
In (a)-(c) three different stable shearons are shown for a chain with free ends. In (d)-(f) periodic boundary conditions are 
used with box sizes i = 13, 13, and 12. Harmonic interaction between the particle is assumed in (a)-(d) and Lennard- Jones 
interaction in (e) and (f). The resulting mean friction forces in (a)-(f) arc (/ k ) 0.263, 0.180, 0.125, 0.123, 0.060, and 0.047. 
The model parameters are a = 1, (3 = 1, 7 = 0.75, A = 0.1, e = 0.01, N = 15, and v = 0.02. [Please note that the figure has a 
reduced quality due to figure size limitations of the condmat database.] 
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FIG. 2. Plot of the particles density g vs position y — y clns and time r for parameters that correspond to periodic stick-slip 
motion. In (a) we focus on the time interval of a single slip event. The slip starts at r = [largest spring force] and persits 
until r = 689 [smallest spring force], (b) and (c) show enlargements of the region where the motion of the top plate starts and 
stops. The model parameters are a = 0.02, j3 = 1, 7 = 0.1, A = 0.1, e = 0.01, iV = 15 [harmonic interaction and free boundary 
conditions], and v = 0.06. [Please note that the figure has a reduced quality due to figure size limitations of the condmat 
database.] 
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FIG. 3. Modifying a shearon by a defect; plotted is the particles density g vs position y — j/ cms and time r. The shallow 
defect is reached by the first particle at r = and left by the last particle at r = 1346. The scattering by the defect increases 
the shearon wave vector, which results in a decrease of the mean friction force from (/ k > ^ 0.263 to (/ k )' 0.178. The 
model parameters are a = 1, (3 = 1, 7 = 0.75, A = 0.1, e — 0.01, iV = 15 [harmonic interaction and free boundary conditions], 
v = 0.02, and h = 0.5. [Please note that the figure has a reduced quality due to figure size limitations of the condmat database.] 
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